Gravitation & Cosmology, Vol. 11 (2005), No. 4 (44), pp. 1-10 
© 2005 Russian Gravitational Society 



CONSTRAINING BOUNCING COSMOLOGY CAUSED BY 
CASIMIR EFFECT 

Wlodzimierz Godlowski^ Marek Szydlowski^ | and Zong-Hong Zhu^ | 

f Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland (WG and MS) 

Marc Kac Complex Systems Research Centre, Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 
(MS) 

X Department of Astronomy, Beijing Normal University, Beijing 100875, China 

We constrain the Friedman Robertson Walker (FRW) model with "radiation-like" contribution to the Friedmann 
equation against the astronomical data. We analyze the observational limitations on a (I + 2:)'' term from supernovae 
type la (SNIa) data, Fanaroff-Riley type lib (FRIIb) radio galaxy (RG) data, baryon oscillation peak and cosmic 
microwave background radiation (CMBR) observations. We argue that it is not possible to determine the energy 
densities of individual components scaling like radiation from a kinematic astronomical test. The bounds for density 
parameter for total radiation-like term can be obtained. We find different interpretations of the presence of scaling 
like radiation term: the FRW universe filled with a massless scalar field in a quantum regime (the Casimir effect), 
the FRW model in a semi-classical approximation of loop quantum gravity, the FRW model in the Randall Sundrum 
scenario with dark radiation or cosmological model with global rotation. In this paper we mainly concentrate on 
the Casimir effect arising from quantum effects of the scalar field. This contribution can describe decaying part 
of cosmological constant. We discuss the back reaction of gravity on Casimir- type force which is a manifestation 
of the vacuum fluctuations of the quantum scalar field at low temperature. It is shown that while the Casimir 
energy gives rise to the accelerating Universe, the cosmological constant term is still required. We argue that a small 
negative contribution of a radiation-like term can reconcile the tension between the observed primordial "^He and 
D abundance. Moreover the presence of such contribution can also remove the disagreement between the Hubble 
parameter Ho values obtained from both SNIa and Wilkinson Microwave Anisotropy Probe (WMAP) satellite data. 
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1. Introduction 

From the recent measurements of distant supernovae 
type la [U [2] we deduce that the universe is in an ac- 
celerating phase of expansion. The effective and sim- 
ple explanation of this current state of the Universe re- 
quires that one third of total energy in the universe is 
non-relativistic dust (dark matter) and two third is a 
constituent of negative pressure (dark energy). While 
cosmological constant remains the simplest explanation 
of observations of distant supernovae, there appeared 
serious problem in this context like: why is the vacuum 
energy so much smaller than we except from effective 
quantum field calculations? 3 . To better understand 
how vacuum energy contributes to the cosmological con- 
stant, we consider the Casimir effect which is a purely 
quantum field theory phenomenon. We discuss the back 
reaction of gravity on Casimir-type force which is a man- 
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ifestation of the vacuum fiuctuations of the quantum 
scalar field at finite temperature. It is shown that while 
the Casimir energy gives rise to accelerating Universe, 
the cosmological constant is still required. The FRW 
model with a Casimir type force contains a term which 
scales like negative radiation ((— )(l-|-z)''). 

There are different interpretations of the presence 
of such a term: cosmological model with global rota- 
tion, Friedmann-Robertson- Walker (FRW) model in the 
Randall Sundrum scenario with dark radiation, FRW 
universe filled with a massless scalar field in a quan- 
tum regime (Casimir effect), or FRW model in a semi- 
classical approximation of loop quantum gravity. To 
constrain the "radiation-like" contribution to the Fried- 
mann equation, we use a variety of astronomical obser- 
vations, such as SNIa data [1, 2 , FRIIb RG data [J], 
baryon oscillation peak and CMBR observations. Al- 
though the obtained bounds on total density parameters 
strongly limit the presence of this term in the Friedmann 
equation, this does not mean that it is not present, as 
it is impossible to determine the energy density of sepa- 
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rate components. We note that the CMBR and big-bang 
nucleosynthesis (BBN) offer stringent conditions on this 
term, which can be regarded as an estabhshed upper 
Hmit on any individual components of negative energy 
density, and therefore on the Casimir effect, global rota- 
tion, discreteness of space following loop quantum grav- 
ity or brane dark radiation. The plan of the paper is 
as follows. In the next section, we summarize differ- 
ent interpretations to the possibility radiation like term 
in Friedmann equation. We provide observational con- 
straints on the radiation like term in Section III. Finally 
we conclude our work in Section IV. 

2. Different interpretations of the 
presence of negative radiation like 
term in the H'^{z) relation. 

2.1. Casimir effect 

Ishak [5D] distinguished old and new cosmological prob- 
lems. While the old problem is related to the order of 
magnitude of cosmological constant the new problem is 
related to this magnitude being of the same order as the 
matter density during the present epoch. He pointed out 
the relevance of experiments focusing on the Casimir ef- 
fect in gravitational and cosmological context [20]. He 
argued that while the geometrical cosmological constant 
has no quantum properties, the vacuum energy has both 
gravitational and quantum properties. In this context 
the Casimir effect [S] seems to be relevant because it can 
tell us how vacuum energy contributes to the cosmolog- 
ical constant [HI [Tj . 

The Casimir effect can be derived from the gen- 
eral principles of quantum theory of electromagnetic 
field [S] for two uncharged, perfectly conducting plates 
in vacuum. They should attract each other with the 
force F cx , where d is distance between the plates. 
This prediction was verified experimentally and yielded, 
for electric bodies, reasonable agreement to the theory 
[8]. The Casimir effect is a simple observational conse- 
quence of the existence of quantum fluctuations [9j . The 
Casimir force between conducting plates leads to a re- 
pulsive force, like the positive cosmological constant. It 
is worthy to mention that Casimir type of contribution 
arising from the tachyon condensation is possible [10] . 

Moreover, different laboratory experiments were de- 
signed to measure the Casimir effect with increased pre- 
cision and then strengthen the constraints on correc- 
tion to Newtonian gravitational law [11]. Therefore, the 
measurement of the thermal Casimir force is promising 
to obtain stronger constraints on non-Newtonian grav- 
ity. Although the Casimir force is very weak (for typical 
values of d = 0.5/im and area of 1 cm^, the Casimir 
force is ~ 0.2 dyn) it becomes measurable with a high 
degree of precision [12] . 

For a survey of recently obtained results in the 
Casimir energy studies see [5]. The Casimir effect 



can be regarded as a manifestation of the quantum 
fluctuations on the geometry and the topology of the 
system boundaries although the Casimir force can be 
also present with the system with no boundaries and a 
compact topology. Therefore, if our Universe has non- 
trivial topology then every quantum flclds will generate 
a Casimir-type force which many authors study on the 
FRW background space-time. If we consider a massless 
scalar field, conformally coupled to gravity, on the back- 
ground of the static Einstein universe then its Casimir 
energy has been shown to have the form a/a^ with the 
value of a = l/(4807r2) [iS;. While for all such fields a 
is positive because such models obey the strong energy 
condition, the computation of Casimir energy in the 
cosmological context leads to the conclusion that the 
Casimir energy of the scalar fields could drive the infla- 
tion in the flat universe with toroidal topology [M]. It 
has been recently demonstrated that there exists a fam- 
ily of quantum scalar fields which give rise to a repulsive 
Casimir force in a closed universe [15] . They "produce" 
Casimir energy scaling like radiation (a < 0), and vio- 
lating the strong energy condition. However, we must 
remember that all calculations of the Casimir effect are 
performed under the assumption of a quasi adiabatic 
approximation and their generalization to the case of 
the non-static FRW models is extremely difficult. 

The finite temperature quantum effects of massless 
scalar fields on the background of spacetimes of cosmo- 
logical models have been considered by many authors 
for many years in the context of the Kaluza-Klein theo- 
ries and dynamical reduction of extra dimension process. 
The main result of these investigations was a universal 
quantum correction in low and high temperatures. As- 
suming the so called adiabatic approximation it can be 
shown that while at low temperature the Casimir energy 
is always negative and proportional to a/a'^'^'^ where 
d is the number of extra dimensions, and is positive 
and proportional to a/a^+'* at high temperatures [TB] . 
Therefore if d = we can recover results for the stan- 
dard cosmology which indicate that quantum effects in 
cosmological models at finite temperatures have a uni- 
versal asymptotic in both high and low temperatures. 
The Casimir energy at low temperatures scales like ra- 
diation but is still negative. Finally the Casimir effect 
is significant when the topology of the Universe is not 
trivial [U]. Also Casimir type energy can be produced 
from some extra dimensions [HI [19] . Recently the rele- 
vance of the Casimir effect in the context of dark energy 
problem has been pointed out [101 dD [H [H [21] . 

When the field occupies some bounded region of the 
configuration space then its spectrum is discrete and 
general vacuum energy Eq depends on eigenfrequen- 
cies ujn which can be determined from the geometry 
of boundary dM ; Eo{dM) — ^J2^n- While the total 
zero-point energy of the vacuum is infinite in the pres- 
ence of boundaries it is modified and Eo — Eo{dA4) — 
£'o(0) In the case considered by Casimir, when electro- 
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magnetic field is confined between two parallel conduct- 
ing plates the contribution of unbounded Minkowski 
space —Eo{0) should be subtracted. The generally to 
obtain the physical value of the vacuum energy we must 
remember that in quantum field theory (i.e., in the case 
of infinite number of degrees of freedom) the observable 
quantity is not zero point energy it self, but only its 
excess, caused by boundaries [S]. 

Recently, the idea that the Casimir effect is respon- 
sible for UV ultraviolet cut-off that renders the total 
vacuum energy finite was rule out ,25j. Therefore the 
Casimir effect cannot be treated as a natural cut off 
leading to the observed cosmological constant value. 

The cosmological significance of the Casimir effect 
was pointed out in many contexts since the pioneer- 
ing paper of Zeldovich and Starobinski jT4]. Also, the 
Casimir effect has been used as an effective mechanism 
of compactification of extra dimensions in the Kaluza- 
Klein cosmological models — so called dynamical reduc- 
tion mechanism ]W[ [55] . 

It was demonstrated that nontrivial topology of a 
physical space as well as an internal space can lead to 
compactification of an extra dimension [27ll28j . (for con- 
temporary context see also [18]. In these investigations 
the quantum field theory at finite temperature is used 
in calculation of quantum effects and back-reaction aris- 
ing from massless scalar fields. It is assumed that the 
metric of background space time (FRW space time) is 
static and conditions for thermodynamical equilibrium 
of matter fields are satisfied due to their interactions 
with the thermal bath. As a result we obtain univer- 
sal approximation of quantum distribution function for 
massless scalar bosons at high and low temperatures. 
The assumption of the quasi-static approximation (char- 
acteristic time of quantum process is much smaller than 
the characteristic time of cosmic evolution) enables us 
to determine the thermodynamical characteristic. As a 
result we obtain at low temperatures: p{M^) — (— )/3, 
p = {—) Jd+4 , p'{M^) — j^p, where p and p' are pres- 
sure on physical {M^) and internal (M^) spaces. D is 
the dimension of space with additional dimensions, a is 
the scale factor. 

Because in this case the energy momentum tensor is 
traceless (massless scalar field) we define the energy mo- 
mentum tensor as = diag(p, — p, — p, — p, — p', — p', —p') 
In the standard space of topology TZ x M.^ , where 
is homogeneous and isotropic space with Robertson- 
Walker metric we know that quantum effects of scalar 
field are equivalent to the effect of fluid with pressure 
p = ^p and p = (— )-^ . It is a universal approximation 
effect of quantum fields originating from massless scalar 
fields at low temperatures (Casimir effect). Analogical 
result was recently obtained by Herdeiro and Sampaio 
[15| . The back-reaction problem is considered by taking 
the semi-classical equation with A 

G^. + Kg^. = T-^"" + < T^u > ■ (1) 



Considering R-W symmetry and matter content in the 
form of perfect fluid with energy density p and pressure 
p we obtain 

a' + k^^, (2) 

d = -i(p + 3p)a, (3) 

where dot means differentiation with respect the cosmo- 
logical time and p is the effective energy density 

PcS = A+^, (4) 

where a is the constant of the Casimir force scaling like 
radiation, positive for conformally coupled scalar held 
in the cosmological context as it was verified long time 
ago by Zeldovich and Starobinski [T3] . Casimir energy is 
negative in a flat universe with toroidal topology. In this 
case ( a < ) it is produced accelerating phase of expan- 
sion of the universe because the strong energy condition 
PcS + 3pcff > is violated. 

Peff-A-- + — , p^, = -A-— (5) 

For our aims it is important that the effects of 
Casimir energy with a negative value of a which scales 
like radiation can contribute into the H'^{z) relation — 
crucial for any kinematic test. It is also interesting that 
the same type of contribution to the effective energy 
density can be produced by loop quantum theory ef- 
fects in semi-classical quantum cosmology 29, ,30] |3I] • 
These effects give rise to evolutional scenario in which 
the initial singularity is replaced by a bounce. 

2.2. Other interpretations of the Casimir-type 
term 

There are many different interpretations of the term in 
the Friedmann equation which diminishes with the cos- 
mic scale factor like a^^ . The first interpretation comes 
from the generalized Friedmann equation on the brane 
in the Randall and Sundrum scenario [35] In the 
brane world scenario our universe is some sub-manifold 
which is embedded in a higher-dimensional space-time 
called bulk spaces. While the physical matter fields are 
confined to this sub-manifold called brane, the grav- 
ity can reside in the higher dimensions. This brane 
paradigm was first proposed by Arkani-Hamed et al. 
[34] as a means to reconcile the hierarchy problem be- 
tween the weak scale and the new Planck scale A^pi. 
Randal and Sundrum [32] solved the analogous problem 
between the weak scale and the size of extra dimension 
by introducing non-compact extra dimensions. In their 
model our universe is represented by a three brane em- 
bedded in a 5 dimensional anti de Sitter space. The 
cosmological evolution of such brane universes was ex- 
tensively investigated by several authors (see for exam- 
ple [35] ). This way, the Einstein equations restricted 
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to the brane reduce to some generalization of the FRW 
equation. Two additional terms contribute into the 
relation \6& . The term arises from the imposition of 
a junction condition for the scale factor on the surface 
on the brane. This term decays rapidly as a^^ for dust 
or as for the radiation dominated early universe. 
This term should be significant only in the very early 
universe [5S1I57] . 

The second term is of considerable interest for us 
because it scales like radiation with a negative constant 
a . Hence it is called dark radiation. This term arises 
from the non-vanishing electric part of the five dimen- 
sional Weyl tensor. Mathematically both negative and 
positive values of a are possible. 

Dark radiation strongly effects both BBN and CMBR. 
It was demonstrated by Ichiki et al. 38J that BBN lim- 
its the possible contribution from dark radiation just 
before e^e~ annihilation epoch. They gave limits on 
the possible contribution of dark radiation as —1.23 < 
Pdrlp-y < 0.11 from BBN and -0.41 < pdr/p-y < 0.105 
at the 95% confidence level from CMBR measurements. 
Let us note that small negative contribution of dark 
radiation can also reconcile the tension between the 
observed ^He and D abundances [39] . 

Another interpretation of the presence of the neg- 
ative radiation like term is rotation in the Newtonian 
cosmology. When we consider Newtonian cosmology fol- 
lowing Senovilla et al. [40] then we can define, homoge- 
neous Newtonian cosmology as p and p having no spa- 
tial dependence i.e p = p(t) and p — p{t) while we as- 
sume that the velocity vector fields depends linearly on 
the spatial coordinates. In such a case we obtain equa- 
tion which represents shear-free Newtonian cosmologies 
with expansion and rotation which is well known as the 
Heckmann-Schiicking model [41] . 



where C is an arbitrary constant. We interpret it in 
terms of curvature constant although in the Newtonian 
spacetime the curvature is zero. For our aims it is im- 
portant that the effect of rotation produce negative term 
scaling like (1 -I- z)^ in the Newtonian analogue of the 
Friedmann equation. 

In the Newtonian cosmology in contrast to general 
relativity effect of rotation are not necessary related to 
non- vanishing shear. The homogeneous universe with 
non- vanishing shear basing on general relativity may ex- 
pand and rotate relative to local gyroscopes. The prob- 
lem of relation between the rotation of the universe and 
origin of the rotation of galaxies was investigated in [42] 
and [43] HH [45]. Also, the role of rotation of objects in 
the Universe, their significance and astronomical mea- 
surements was recently addressed by [46l [47] . 

To compare the results of analyses for general rela- 
tivistic model with results obtained further in the paper 
we formally consider Qk,o 7^ 0, although the satisfac- 
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Figure 1: The 68.3% and 95.4% confidence levels (obtained 
from combined analysis of SN-^RG-hSDSS-hCMBR) on the 
(rim.o, f^A.o ) plane. 
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Figure 2: The location of the first peak h as a function of 
—Qdr,o- Note that h — 220 for h — 0.73 favour ^ldr,o — 0, 
while for h = 0.65 and h — 0.67.6 lead to fldrfi 7^ 0. 



tory interpretation of curvature density parameter ilkfi 
can be found in general relativity (see also [37]). 

3. Observational constraints on the 
FRW model parameter with the 
radiation-like term 

Cosmological model are usually tested against observa- 
tions. One of the most popular test is based on the 
luminosity distance c?l of the supernovae la as a func- 
tion of redshift [H] . However for the distant SNIa, one 
can not directly observe their luminosity distances cIl 
but their apparent magnitudes m and redshifts z . The 
absolute magnitude A4 of the supernovae is related to 
its absolute luminosity L. Since we could obtain the 
following relation between distance modulus p, the lu- 
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minosity distance , the observed magnitude m and 
the absolute magnitude M : 



M 



5 logio ' 



-Blog^Q Dl+M (7) 

~51ogioi7o + 25. We 
could compute the luminosity distance of a supernova 
as the function of rcdshift: 



where Dl — HadL and Ai 



dL{z) = (1+z) 



1 



Ia\J\^kfl\ 



dz' 




(81 



where 



— ) = a„.o(l + z)^ + f^fe,o(l + z)^ + f^,,o(l + z 

^0 



— -ffi and J-{x) — (sinh(a;), sin(a;)) for k 
,k>0, respectively. We assumed il^.o — ^ 



< 



7,0 ' 



^k,0 ~ 

0, fc = 

r2^,o = 2.48/i-^ X 10"^ + 1.7h-'^ X 10"'^ ~ 0.0001 [36 

Daly and Djorgovski [49] suggested to include in the 
analysis not only supernovae but also radio galaxies (see 
also [ini [SU [52]). In such a case, it is useful to use the 
coordinate distance defined as 



viz) 



c(l + z) ■ 



(10) 



Daly and Djorgovski [4] have compiled a sample com- 
prising the data on y{z) for 157 SNIa in the Riess et al. 
Gold dataset [I] and 20 FRIIb radio galaxies. In our 
data sets we also include 115 SNIa compiled by Astier 
et al. [2]. 

It is clear from eq. © and pU|) that the coordinate 
distance does not depend on the value of iJp • Unfor- 
tunately we do not know coordinate distance y{z) for 
supernovae. This distance must be computed from the 
luminosity distance (or the distance modulus /i ) and for 
such a computation a knowledge of the value of Hq is 
required. For both supernovae sample we choose the 
values of which were used in the original papers. 
We used the distance modulus presented in Ref. [2 [5] 
for the calculation of the coordinate distance. For each 
sample we choose the values of Hq appropriate to the 
data sets. For Riess et al.'s Gold sample we fit the value 
oi h — 0.646 as the best fitted value and this value is 
used for calculation of coordinate distance for SNIa be- 
longing to this sample. In turn the value h — 0.70 was 
assumed in the calculations of the coordinate distance 
for SNIa belonging to Astier et al.'s sample, because the 
distance moduli /i presented in Ref. [2j Tab. 8] was cal- 
culated with such an arbitrary value oi h = 0.70. The 
error of the coordinate distance can be computed as 



10- 



c (1 + z) 105 



X a\Ho) 



'i/oinioy 2 



'^'{^^^) \ (11) 



where cri{fii) denotes the statistical error of distance 
modulus determination (note that for Astier et al.'s 
sample the intrinsic dispersion was also included) and 
iT(-ffo) = 0.8 km/s Mpc denotes error in Hq measure- 
ments. 

We included to our constraints obtain from extra- 
galactic analysis: the baryon oscillation peaks (BOP) 
detected in the Sloan Digital Sky Survey (SDSS) [53]. 
They found that value of A 



1 



£■(21)3 \ziy/\nkfi\ 



fe,0| 



dz 



(12) 



(where E{z) = H{z)/Hq and zi — 0.35) is equal A = 
0.469 ± 0.017. The quoted uncertainty corresponds to 
one standard deviation, where a Gaussian probability 
distribution has been assumed. 

Another constraint which we also include in our anal- 
ysis is the so called the (CMBR) "shift parameter" 



dz 



R = \J^iafiy{ziss) 




E{z) 



where Rq = 1.716 ± 0.062 

In our combined analysis, we can obtain a best fit 
model by minimizing the pseudo-x^ merit function [55j 

2 _ 2 I 2 I 2 

X — XSN+RG + XSDSS + XCMBR 



E 



obs _ til 



^mod _ 0.469 
0.017 



1.716 



0.062 



(14) 



(13) 



where A'^°'^ and denote the values of A and R 

obtained for a particular set of the model parameter. 
For Astier et al. SNIa [2] sample additional error in 
z measurements were taken into account. Here cri(%ji) 
denotes the statistical error (including error in z mea- 
surements) of the coordinate distance determination. 

Constraints for the cosmological parameters, can be 
obtain by minimizing the following likelihood function 
C cx exp(— x^/2) . One should note that when we are 
interested in constraining a particular model parameter, 
the likelihood function marginalized over the remaining 
parameters of the model should be considered [55] . This 
method was used in the case of negative (1 + z)'^ type 
contribution only in the paper [52 . In the present paper 
we avoid from constrain that a priori the (1 + z)^ term 
must be negative. 

Our results are presented in Table [TJ Table [2| and 
Fig. [H Table [T] refers to the minimum method, 
whereas Table [31 shows the results from the marginal- 
ized likelihood analysis. 

From our combined analysis (SN+RG+SDSS+CMBR) 
we obtain as the best fit a flat (or nearly flat universe) 
with r^ni^o — 0-3, and JIa.o — 0.7. For the dark ra- 
diation term, we obtain the stringent bound f^Totr.o = 
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^r,o + ^dr.o e (-0.00017,0.00857) at the 95% confi- 
dence level. It lead for limit on dark radiation ildr.o G 
(-0.00027,0.0847) (f^Totr.o e (-0.00017,0.0857). This 
results mean that the positive value of dark radiation 
term is preferred {^dr,o > 0), however small negative 
contribution of dark radiation is also available. One 
should note that when we used only SN and RG data 
only, we obtain value Qrii,o close to zero what seems 
to be unrealistic. It is the reason that we repeat our 
analysis with prior ilm.o = 0.3 [56]. In such a case, 
we again obtain flat (or nearly flat universe) as a best 
fit. Contribution of dark energy is small and positive, 
but also in this case small negative contribution of dark 
radiation is available. From combined analysis we ob- 
tain Hmit on dark radiation f7dr,o e (-0.00037, 0.00727) 
(f^Totr,o e (-0.00027,0.00737)). Our results shows that 
in the present epoch contribution of the dark radiation, 
if it exist, is small and gives only small corrections to 
the ACDM model in the low redshift. 

We use the Akaike information criteria (AIC) [60] 
and the Bayesian information criteria (BIG) fEV to se- 
lect model parameters providing the preferred fit to 
data. Usually incorporating new parameters increases 
the quality of the fit. The question is, if it increases 
it significantly enough. The information criteria put a 
threshold which must be exceeded in order to assert an 
additional parameter to be important in explanation of 
a phenomenon. The discussion how high this thresh- 
old should be caused the appearing of many different 
criteria. The AIC and BIG (for review see [ST]) are 
most popular and used in everyday statistical practices. 
The question is whether the AIC and BIG information 
criteria of model selection provide sufficient arguments 
for incorporation of new parameters. The usefulness 
of using the information criteria of model selection was 
demonstrated by Liddle [58l and Parkinson et al. [52 Of 
course, in any case, some future observational data, as 
well as theoretical considerations, may give arguments 
in favour of additional parameter. 

The AIC is defined in the following way [50] 

AIC = -21n/: + 2d (15) 

where £ is the maximum likelihood and d is a number 
of the free model parameters. The best model with a pa- 
rameter set providing the preferred fit to the data is that 
minimizes the AIC. It is interesting that the AIC also 
arises from an approximate minimization of the Kulbak- 
Leibner information entropy [62] . 

The BIG introduced by Schwarz [5T] is defined as 

BIG = -21n/: + <ilniV (16) 

where N is the number of data points used in the fit. 
Comparing these criteria, one should note that the AIC 
tends to favour models with large number of parameters 
unlike the BIG, because the BIG penalizes new param- 
eters more strongly. It is the reason that the BIG pro- 
vides a more useful approximation to the full statistical 



analysis in the case of no priors on the set of model pa- 
rameters [59] . It makes this criterion especially suitable 
in the context of cosmological applications. 

Please note that while the AIC is useful in obtaining 
upper limit to the number of parameters which should 
be incorporated to the model, the BIG is more conclu- 
sive. Of course only the relative value between the BIG 
of different models has statistical significance. The dif- 
ference of 2 is treated as a positive evidence (and 6 as a 
strong evidence) against the model with the larger value 
of the BIG [63l [64] . If we do not find any positive evi- 
dence from information criteria the models are treated 
as a identical and eventually additional parameters are 
treated as not significant. The using of the BIG seems 
to be especially suitable whenever the complexity of ref- 
erence does not increase with the size of data set. Liddle 
[58] noted that in cosmology, a new parameter is usu- 
ally a quantity set to zero in a simpler base model and 
if the likelihood function is a continuous function of its 
parameters it will increase as the parameter varies in ei- 
ther the positive or negative direction. The problem of 
classification of the cosmological models on the light of 
information criteria on the base of the astronomical data 
was discussed in our previous papers [65l [66l [Ml [69] • 

Our results are presented in Table [5j It is clear that 
in the light of informative criterion model with dark 
energy do not increase fit significantly. It confirm our 
conclusion that dark energy term, if it exist, is small in 
the present epoch. Using the prior flni,o = 0.3 [56] does 
not change our conclusion. It shows that using such a 
prior in the light of information criteria is realistic. 

Please also note that if riTotr,o = ^i,o + ^dr,o < , 
then we obtain a bouncing scenario [TQl [Til [22] instead 
of a big bang. For fJm.o = 0.3, fldr.o = —0.00027 and 
h = 0.65 bounces {H^ — 0) appear for z ~ 1800. In 
this case, the BBN epoch never occurs and all BBN 
predictions would be lost. 

This results shows that more stronger constraints for 
model parameters is required. To obtain such a con- 
straints on the model parameters, it is useful to use the 
CMBR observations. The hotter and colder spots in 
the CMBR can be interpreted as acoustic oscillations in 
the primeval plasma during the last scattering. In such 
a case the locations of the peaks in the CMBR power 
spectrum are very sensitive to variations in the model 
parameters. Therefore, the location of these peaks can 
also be used for constraining the parameters of cosmo- 
logical models. The acoustic scale £a which gives the 
locations of the peaks is defined as 

r^doc dz' 

where, for the fiat model, equation ^ reduces to 

H{z) = i/o [f^m,0(l + ^)' + r!r.0(l + ^)^ 

+ nd.fi{l + zY + n^fl]''\ (18) 
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where Cg is the speed of sound in plasma. Knowing 
the acoustic scale we can determine the location of the 
m-th peak £m ~ — (I'm) where (pm is the phase 

shift caused by the plasma driving effect. The CMBR 
temperature angular power spectrum provides the lo- 
cations of the first two peaks ii = 220.llQg, £2 = 
546ti!^ IZ3J. Using three years of WMAP data, Spergel 
et al. obtained that the Hubble constant Hq = 73 km/s 
Mpc, the baryonic matter density fib.o = 0.0222/i^^, 
and the matter density ^1^,0 = 0.128/i~^ [74j, which 
are in good agreement with the observation of position 
of the first peak (see Fig. [2]) but lead (assuming the 
ACDM model) to a value rim.o = 0.24. It mean that 
there is disagreement between Hq values obtained from 
SNIa and WMAP. We compute the location of the first 
peak as a function of Odr.o assuming Hq — 65 km/s 
Mpc (ilm.o = 0.3). Separately we repeat our compu- 
tation using the latest Riess et al. result which obtain 
n„,,o = 0.28 (Ho =67. 6 km/s Mpc) ^75^. 

From Fig. [2] it is easy to see that we can obtain 
agreement with the observation of the location of the 
first peak for non-zero values of the parameter Qdi,a 
(Fig. ^ both for n„,^o = 0.28 and f7n,,o = 0.30. Wc 
obtain -1.05 x 10^^ < ridr,o < -0.5 x lO"'^ at the 
95% confidence level for the case rim.o — 0.3 while 
-0.75 X 10-^ < 17dr,o < -0.25 x lO^^ at the 95% 
confidence level for the case r2m,o = 0.28 Please note 
that our limits are stronger than that obtained by Ichiki 
et al. [38], which provides bounds of —7.22 x 10^''' < 
f^dr.o < 0.65 X 10^5 (in the case of the BBN) and 
-2.41 X 10"^ < rtdrfi < 0.62 X 10"^ (in the case of 
the CMBR). 

In all cases the obtained values of fidr.o a-re in agree- 
ment with the result obtained from the combined analy- 
sis because the 2a confidence interval for this parameter 
obtained from this analysis contains the area allowed 
from the CMBR. Most important conclusion is, that 
while the combined analysis allowed the possibility that 
^^dr,o is equal to zero, the CMBR location of the first 
peak seems to exclude this case both for h = 0.65 and 
h = 0.67. 

4. Conclusion 

In the paper we analysed the observational constraints 
on the (l-|-2:)*-type contribution in the Friedmann equa- 
tion. Because it the present paper we mainly concen- 
trate on Casimir effect arising from quantum effects of 
the scalar field the constraints for negative (l-l-z)^-type 
contribution are in our special interest. The analysis of 
SNIa data as well as both SNIa and FRIIb radio galaxies 
(with and without priors going from baryon oscillation 
peaks and CMBR "shift parameter" ) shows that the val- 
ues of statistics are lower for model with radiation 
like term, than for the ACDM model. However, infor- 
mation criteria show that using such a term does not 
increase the quality of the fit significantly. BIC even 



favour A CDM model over our model of Bouncing Cos- 
mology (with dark radiation). However this preference 
is weak. This results show that (1 -I- z)^ term is not 
significant in the present epoch of the Universe. 

We show that there are several interpretations of 
the (1 -I- z)^-type contribution and we discussed dif- 
ferent proposals for the presence of such a term. Un- 
fortunately, it is not possible, with present kinematic 
astronomical test, to determine the energy densities of 
individual components scales like radiation. However 
we show that some stringent bounds on the value of 
this total contribution can be given. Combined analysis 
of SNIa data and FRIIb radio galaxies using baryon 
oscillation peaks and CMBR "shift parameter" give rise 
to a concordance universe model which is almost fiat 
with rim.o — 0.3 . From the above-mentioned com- 
bined analysis, we obtain an constraint for the term 
which scales like radiation ilTotr,o e (-0.00017,0.00857) 
which leads to bounds on the dark radiation term 
f^dr,o e (-0.00027,0.0847) at the 95% confidence level. 
This is a stronger limit than obtained previously by us 
from SNIa data only [47J. Our model with a small con- 
tribution of dark radiation type can also resolve the dis- 
agreement between Hq values obtained from SNIa and 
WMAP. For 17ni,o = 0.3 (ffo = 65 km/s Mpc) we find 
new stringent limits on the negative component scaling 
like radiation from the location of the peak in the CMBR 
power spectrum, —1.05 10^^ < ildr,o < —0.5 10~^ at 
the 95% confidence level. This bound is stronger than 
that obtained from BBN and CMBR by Ichiki et al. 

In this paper we have especially studied the advance 
of initial singularity using back reaction gravity quan- 
tum effect at low temperatures (Casimir effect). Casimir 
force arising from the quantum effect of massless scalar 
field give rise to a (— )(1 -I- z)^ correction whose effect 
depends upon the geometry and nontrivial topology of 
the space. Typically this type of correction is thought 
to be important at the late time of evolution of the uni- 
verse. We have shown that Casimir effect could remove 
initial singularity which would be replaced by bounce. 
However, from the observational limits we obtain that 
bounce does not appear which means that a big bang 
scenario is strongly favoured instead of bounce. Our 
limit {—)^ldr,o < 1-05 10~^ obtain from position of the 
first peak in the power spectrum of CMBR leads to 
^^Totr.o — ^dr.a + ^r,o > 0. This implies that H'^{z) 
is always greater than zero {H^{z) > 0) and bounce 
does not appear which means that a big bang scenario 
is strongly favoured instead of bounce. 
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Table 1: Results of the statistical analysis of the model with radiation like term obtained from best fit. The 
upper section of the table represents the constraint Qk,o = (flat model). 



sample 






^^Totr,0 






SN 




0.01 


0.133 


0.807 


295.9 


SN+RG 




0.11 


0.123 


0.767 


319.5 


SN+RG+SDSS 




0.28 


0.022 


0.698 


320.1 


SN+RG+SDSS+CMBR 




0.30 


0.00054 


0.699 


322.3 


SN 


0.01 


0.00 


0.186 


0.804 


295.9 


SN+RG 


0.08 


0.00 


0.169 


0.751 


319.5 


SN+RG+SDSS 


-0.11 


0.28 


0.049 


0.781 


319.5 


SN+RG+SDSS+CMBR 


0.03 


0.29 


0.00270 


0.657 


320.9 



Table 2: Results of the statistical analysis of the model with radiation-like term. The values of the model parameters 

arc obtained from marginalized likelihood analysis. Wc present maximum likelihood value with 68.3% confidence 
ranges. The upper section of the table represents the constraint flkfi = (flat model). 



sample 



k,0 



SN 

SN+RG 

SN+RG+SDSS 

SN+RG+SDSS+CMBR 



SN 

SN+RG 

SN+RG+SDSS 

SN+RG+SDSS+CMBR 



^•^"-0.28 

-0.19+0-24 



^•"^-0.00 
-,+0.16 
-0.10 

h0.02 
-0.01 



0.10] 



0.30 



i^Totr.O 

0.045 
-0.263 
0.057 
0.109 
0.015 
0.014 



0.164 
0.128 
0.022 



0.27 
+0.15 
-0.16 

06+0 °^ 



-0.11 



hO.bV 
-0.00 
hO.55 
-0.00 
h0.02 
-0.02 

n 28+°-°2 



0.00056"JS5 



0.00] 
0.00^ 
0.28^ 



0.091 



h0.0b3 
-0.201 

U.U01_Q 
nm-^A+O -00267 
U.UUJJ4t_Q QQ226 



f^A,0 

n 70 + 0.03 

*J-'^-0.08 
n 77+0.03 
'-0.08 

"•"^-0.01 

0.811°:- 

0.781°:- 
0-64tro^ 



Table 3: Results of the statistical analysis of the model with the radiation like term obtained from best fit with 
the assumption Om,o = 0.3. The upper section of the table represents the constraint flk,o = (flat model). 



sample 




f^m.O 


f^Totr,0 


^A,0 


x' 


SN 




0.30 


0.009 


0.691 


297.3 


SN+RG 




0.30 


0.010 


0.690 


320.2 


SN+RG+SDSS 




0.30 


0.012 


0.688 


321.0 


SN+RG+SDSS+CMBR 




0.30 


0.00054 


0.699 


322.3 


SN 


-0.19 


0.30 


0.056 


0.834 


295.9 


SN+RG 


-0.13 


0.30 


0.042 


0.788 


319.5 


SN+RG+SDSS 


-0.09 


0.30 


0.034 


0.756 


320.7 


SN+RG+SDSS+CMBR 


0.03 


0.30 


0.00182 


0.668 


321.4 



Table 4: Results of the statistical analysis of the model with the radiation like term. The values of the model 

parameters are obtained from the marginalized likelihood analysis with the assumption f^m.o = 0.3. We present 
maximum likelihood value with 68.3% confidence ranges. The upper section of the table represents the constraint 
^fe,o = (flat model). 



sample 


^k,0 f^m.O 


f^Totr.O 


^^A,0 


SN 

SN+RG 


0.30 
0.30 


009+° '^^^ 

u.uuy_Q 010 
010+° °" 


69+° °^ 
0.69t°:°^ 


SN+RG+SDSS 


0.30 


012+°-°i° 


0.68t-°? 


SN+RG+SDSS+CMBR 


0.30 


n ooo^4+°°°i62 

U.UUUO+_o. 00070 


0.69t<°:°} 



SN 


-0 ig+o-^'' 

"■-'-^-0.17 


0.30 


n's7+"-°'^^ 


0.87t°:l^ 


SN+RG 


-0 1S+°15 


0.30 


042+° °''° 


0.82t°:- 


SN+RG+SDSS 


-0 09+°-^^ 

"•'^^-0.15 


0.30 


034+° 039 


o.6ii°:- 


SN+RG+SDSS+CMBR 


04+° °^ 

'-'•"^-0.03 


0.30 


00217+° °°237 

U.UU/i r_o. 00173 


0.65;°:°^ 
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Tabk^ 5: The' vahu's of tlic^ AIC and BIC for XXiq ACDM modc^l and Boundng Cosmology model (with dark radiation) 
without (Uid with \moY i^ni.o = Th(> upper section of \ rt^prt^sents the constraint Q-kfi = (flat model). 



ACDM BC BC(Om,o = 0-3) 



sample 


AIC BIC 


AIC BIC 


AIC BIC 


SN 
SN+RG 
SN+RG+SDSS 
SN+RG+SDSS+GMBR 


299.5 303.1 
322.4 326.1 

324.4 328.1 

324.5 328.2 


299.9 307.1 
323.5 330.9 

324.1 331.5 
326.3 333.7 


299.3 302.9 

322.2 325.9 
323.0 326.7 

324.3 328.0 


SN 
SN+RG 
SN+RG+SDSS 
SN+RG+SDSS+CMBR 


300.0 307.2 
323.5 330.9 

325.1 332.5 
326.5 333.9 


301.9 312.7 

325.5 336.5 
325.5 336.5 
326.9 337.9 


299.9 307.1 

323.5 330.9 
324.7 332.1 
325.4 332.8 



